Random walks on the Apollonian network with a single trap 
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Abstract. - Explicit determination of the mean first-passage time (MFPT) for trapping prob- 
lem on complex media is a theoretical challenge. In this paper, we study random walks on the 
Apollonian network with a trap fixed at a given hub node (i.e. node with the highest degree), 
which are simultaneously scale-free and small-world. We obtain the precise analytic expression for 
the MFPT that is confirmed by direct numerical calculations. In the large system size limit, the 
MFPT approximately grows as a power-law function of the number of nodes, with the exponent 
much less than 1, which is significantly different from the scaling for some regular networks or 
fractals, such as regular lattices, Sierpinski fractals, T-graph, and complete graphs. The Apol- 
lonian network is the most efficient configuration for transport by diffusion among all previously 
studied structure. 



Introduction. — Trapping is an integral major theme 
of random walks (diffusion) [1-3], which is relevant to a 
wide range of apphcations and has led to a growing num- 
ber of theoretical and practical investigation over the past 
decades [4-9]. The trapping problem, first introduced 
in [10], is in fact a random-walk issue, where a trap is posi- 
tioned at a given location, which absorbs all particles vis- 
iting it. The primarily interesting quantity closely related 
to trapping problem is the average trapping time, also re- 
ferred to as the mean first-passage time (MFPT) , which is 
useful in the study of transport-limited reactions [11, 12], 
target search [13, 14] and other physical problems. 

An important question in the study of trapping is how 
the MFPT scales with the size of the system. There are 
some well-known results providing answers to the corre- 
sponding c^ucstions in the cases of some graphs with simple 
topology, including regular lattices [10], Sierpinski frac- 
tals [15, 16], T-fractal [17], and so on. However, these 
graphs are not suitable to describe real systems [18] en- 
countered in everyday experience, most of which are scale- 
free [19] and small-world [20] that have been shown to 
influence profoundly various dynamical processes running 
on networks [21,22]. Thus, it is natural and interesting 
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to explore the trapping problem on networks with gen- 
eral structure embedded in real life. Although a lot of 
activities have been devoted to studying random walks on 
complex networks [23-27], work about trapping problem 
on scale- free small- world graphs is much less [28] . 

In the paper, we investigate the trapping problem on 
the ApoUoian network [29, 30] with scale- free and small- 
world properties. We focus on a specific aspect of random 
walks in the presence of a single trap situated at a given 
node with the largest degree (hub node). We obtain an 
exact analytical solution for the MFPT and the depen- 
dence of this primary quantity on the system size. We 
show that the ApoUoian network is a preferred architec- 
ture that minimizes the increase of MFPT with network 
size, compared with regular networks and fractals. 

Introduction to the Apollonian network. We 

first introduce the Apollonian packing [31], from which the 
Apollonian network is derived. There are two commonly 
used Apollonian packings that differ mainly in initial con- 
figurations. The first packing is constructed by starting 
with three mutually touching disks, the interstice of which 
is a curvilinear triangle. In the first generation a disk is 
inscribed, touching all the sides of this curvilinear trian- 
gle. For subsequent generations we indefinitely repeat the 
packing process for all the new curvilinear triangles. In 
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Fig. 1: (Left panel) The first three generations of the first Apol- 
lonian packing of disks. (Right panel) An Apollonian network 
corresponding to the packing shown in the left panel. 

the limit of infinite generations, an Apollonian packing is 
obtained. The left panel of figure 1 displays the first three 
generations of this Apollonian packing. 

The other frequently used Apollonian packing has a ini- 
tial configuration with three mutually touching disks in- 
scribed inside a circular space. The interstices of the initial 
disks and circle are four curvilinear triangles to be filled. 
Then in each subsequent generation, we add one disk to 
each interstice (curvilinear triangle), so that the added 
disk touches all the three sides of the corresponding curvi- 
linear triangle. The first several processes are shown in 
figure 2. 

The above two Apollonian packings can be easily 
mapped to networks, usually called as Apollonian net- 
works [29, 30]. The translation from Apollonian pack- 
ing construction to Apollonian network generation is quite 
straightforward. Each node (vertex, site) corresponds to 
a disk, and two nodes are linked to each other if their cor- 
responding disks are tangent. Note that for the second 
Apollonian packing, the initial circle also corresponds to a 
node. The right panel of figure 1 shows an example of the 
network. Since the resulting networks associated with the 
two packings have similar structural features, except for 
the initial disks, in this paper we will focus on the network 
based on the second Apollonian packing, which is conve- 
nient for analytically deriving the network properties. 

According to the foregoing mapping, one can introduce 
a general algorithm [32,33] to create the Apollonian net- 
work, denoted by Ag after g generation evolutions. For 
5( = 0, Aq is a tetrahedron with four faces or triangles, see 
figure 3. For g > 1, Ag is obtained from Ag_i. For each of 
the existing triangles of Ag_i that is created at generation 
i — 1, a new node is added and connected to all the three 
nodes of this triangle. Figure 4 illustrates the construc- 
tion process for the first two generations of the initial four 
faces as shown in figure 3. 

According to the network construction, one can see [30, 
33] that at the gith {gi > 1) generation, the number of 
newly introduced nodes is Ly{gi) = 4 x 3^'"^. From this 
result, we can easily compute the total number of nodes 
Ng in network Ag (hereafter called the network order): 

a 

Ng=Y,Lv[g^) = 2y.Z3 + 2. (1) 




Fig. 2: The first several generations of the second Apollonian 
packing of disks within a circle. 



For convenient use in the next section, we distinguish 
different nodes of Ag by labeling them sequentially as 
follows. The four nodes (i.e., nodes B, C, D, and E 
in figure 3) created at initial generation are labeled as 
1, 2, 3, and 4, respectively. Then, in each new gener- 
ation, we only label the new nodes added at this gen- 
eration, while keep the labels of pre-existing nodes un- 
changed. That is to say, we label new nodes progressively 
as M + 1,M + 2, . . . ,M + AM, in which M is the total 
number of the pre-existing nodes and AM is the number 
of newly-added nodes. Eventually, every node is labeled 
by a unique integer, all nodes in Ag are labeled from 1 to 
Ng, see figure 4. 

Let ki(g) be the degree of a node i at generation g, 
which entered the network at the generation of gi (gi > 0). 
Then [30,33] 

hig) ^3x2^-3'. (2) 

From equation (2), one can easily see that at each step the 
degree of a node doubles, i.e., 

h{g) = 2h{g-l). (3) 

The Apollonian network presents the typical character- 
istics of real-life networks [29,30,32-34]. It has a power- 
law degree distribution P{k) ~ k~'^ with the exponent 
7=1-1- 1^ . Its average path length, defined as the mean 
of shortest distance between all pairs of nodes, increases 
logarithmically with network order [34] . In the large net- 
work order limit, the average clustering coefficient tends 
to 0.8284. Thus, the Apollonian network exhibits small- 
world effect [20]. In addition, the network is disassor- 
tative [30]. In view of its structural properties similar to 
those of real networks and its intrinsic interest, Apollonian 
network is a good substrate network for studying critical- 
ity phenomena and dynamical processes [35-43]. In what 
follows we will investigate the MFPT for random walks 
with an immobile trap on the Apollonian network. 

Formulation of the problem. In this section 
we formulate the problem of a simple random walk of a 
particle on Apollonian network Ag in the presence of a 
trap or a perfect absorber positioned on a given node. 
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Fig. 3: The initial construction of Apollonian network. 

To this end, we first speeify Ag by its adjaeeney ma- 
trix Ag (hereafter we neglect the subscript) of order 
Ng X Ng, which completely describes network Ag. The 
element ay of A is defined as follows: = 1 if the 
pair of node i and node j is connected by a link (edge, 
bond), otherwise = 0. The degree ki{g) of node i 
is ki (g) = X] j^i '^ij ' and the diagonal degree matrix Z 
is given by Z = diag(fci(g), k2{g), ■ ■ ■ , hig), ■ ■ ■ , kN,{g))- 
Then, the normalized Laplacian matrix of Ag is given by 
L = I — Z~^A and its entry is lij = 1 — jt^i where I is 
an identity matrix. 

Before proceeding further, let us introduce the so-called 
simple random walk [1-3] on network Ag. At each time 
step (taken to be unity), the walker moves from its current 
location to any of its nearest neighbors with equal proba- 
bility. According to this rule, at time t, a walker at a node 
i will hop to one of its ki{g) neighbors, say u, with the 
transition probability aiu/ki{g). Suppose that the walker 
starts from node i at t ~ 0, then the jumping probabil- 
ity Pij of going from i to j at time t is governed by the 
following master equation: 

w — l 

We next focus a random walk on Ag with a trap. We 
locate the trap at node 1, represented as it- It should be 
pointed out that thanks to the symmetry, the trap can be 
also situated at node 2, 3, or 4, which has no influence 
on MFPT. The particular selection for the trap location 
allows one to calculate analytically the MFPT, which will 
be addressed in detail in the next section. At each time 
step, the walker, starting from any node except the trap 
ir, jumps to any of its nearest neighbors with equal prob- 
ability. 

Let Ti be mean transmit time (first-passage time, trap- 
ping time, or the mean time to absorption) for a walker, 
starting from node i, to first arrive at the trap ix- Thus, 
Ti^ = 0. Then, the set of these interesting quantities fol- 
lows the recurrence equation [44] 

where i ^ it- Equation (5) may be also rewritten in 
matrix notation as 

T = A-ie, (6) 
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Fig. 4: The Apollonian network of generation 2 and the labeling 
of its nodes. 

where T (T2, T3, • • • , T^v)^ and e = (1,1,---,1)^ are 
two {Ng — l)-dimensional vectors, and is the fun- 

damental matrix of the Markov chain representing such 
unbiased random walk. In fact, A is a sub-matrix of the 
normalized discrete Laplacian matrix L whose first row 
and column, corresponding to the trap node, have been 
removed. 

Then, the mean first-passage time, or the average of the 
mean time to absorption, (T)g, which is the average of Ti 
over all nodes distributed uniformly over nodes in Ag other 
than the trap, is given by 

^ i=2 y 1=2 j=2 

Equation (7) can be easily explained by looking at the 
random walk from the perspective of a Markov chain. Ac- 
tually, the entry {A~^)ij of the fundamental matrix A~^ 
of the Markov process expresses the average number of 
times that a walker starting at node i will be at node j, 
and the row sum J2j^=2 is exactly Ti, the total 

times a particle starting at node i will traverse all other 
nodes before being absorbed by the trap. 

Equation (7) shows that the problem of determining 
(T)g is reduced to finding the sum of all elements of the 
bmdamental matrix A^^ of order {Ng — 1) x {Ng — 1), 
which can be obtained by utilizing a standard software 
package, Mathematica 5.0. Irrespective of the seemingly 
compact expression of equation (7), since Ng increases ex- 
ponentially with g, for large it becomes difficult to ob- 
tain (T)g through direct calculation from this equation, 
because of the limitations of time and computer memory. 
Therefore, one can compute directly the MFPT only for 
the first generations, see table 1. However, the particular 
construction of the Apollonian network allows one to cal- 
culate analytically MFPT to obtain a rigorous solution. 
Details of derivation will be provided below. 

Exact solution for mean first-passage time. 

Before deriving the general formula for MFPT, (T)g, 
we first establish the scaling relation dominating the evo- 
lution of Tf with generation g, where T[ is the trapping 
time for a walk originating at node i on the gth generation 
of Apollonian network. 
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Table 2: Mean time to absorption Tf for a random walker starting from node i on the Apollonian network for various g. Notice 
that owing to the obvious symmetry, nodes in a parenthesis are equivalent, since they have the same trapping time. All the 
values are calculated straightforwardly from equation (7). 



(2,3,4) 



(5,6,7) 



(8) 



(9,10,12,13,15,16) 



(11,14,17) 



(18,19,20) 



3 

27/5 
243/25 
2187/125 
19683/625 
177147/3125 



23/5 
207/25 
1863/125 
16767/625 
150903/3125 



32/5 
288/25 
2592/125 
23328/625 
209952/3125 



1594323/15625 1358127/15625 1889568/15625 



7 256/25 283/25 

63/5 2304/125 2547/125 

567/25 20736/625 22923/625 

5103/125 186624/3125 206307/3125 

45927/625 1679616/15625 1856763/15625 



Table 1: The MFPT obtained by direct calculation from equa- 
tion (7) by using a standard software package, Mathematica 
5.0. Since for large networks, the computation of the MFPT 
from equation (7) is prohibitively time and memory consuming, 
we calculate the MFPT for the first several generations. 



9 N, 

4 9/3 

1 8 182/(5 X 7) 

2 20 861/(5 x 19) 

3 56 109854/(125 x 55) 

4 164 2895129/(625 x 163) 

5 488 77327622/(3125 x 487) 
6 1460 415448109/(3125x 1459) 



Evolution scaling for trapping time. We begin by 
recording the numerical values of Tf . Obviously, for all 
gr > 0, rf = 0; for g = 0, it is a trivial case, we have 



T2 = T3 = T4 = 3. For .9 > 1, the values of Tf can be 
obtained straightforwardly via equation (7). Table 2 lists 
the numerical values of Tf for some nodes up to g = 6. 
The numerical values listed in table 2 show that for a 
given node i we have Tf^^ — ^Tf. That is to say, upon 
growth of the Apollonian network from g to generation 
g + I, the mean time to first reach the trap increases by 
a factor |. This is a basic character of random walks on 
the Apollonian network, which can be established from 
the arguments below [1,25]. 

Consider an arbitrary node i in the Apollonian network 
Ag after g generation evolution. From equation (2), we 
know that upon growth of the Apollonian network to gen- 
eration 5 + 1, the degree ki of node i doubles. Thus, at 
generation g + 1, node i has ki old neighbors (born at gen- 
eration g or earlier) and ki new neighbors (generated at 
generation g + 1). Let the mean transmit time for going 
from node i to any of its ki old neighbors be X; and let 
the mean transmit time for going from any of its ki new 
neighbors to one of the ki old neighbors be Y. Then we 
can establish the following underlying backward equations 




Fig. 5: Growth of trapping time in going from Ag to Ag+i. 
Node i G Ag has ki neighbor nodes in generation g (■) and ki 
new neighbor nodes in generation g + 1 (□). A new node has a 
degree of 3, and is linked to i and its two old neighbor nodes. 

(see figure 5) 

X = ^ + ^il + Y), y^^ + i(i + x), (8) 

which leads to X = |. That is to say, the passage time 
from any node i (i € Ag) to any node j (j S Ag) increases 
by a factor of | , upon the network growth from generation 
g to generation g + 1. Thus, we have Tf^^ ~ $ Tf , which 
will be useful for deriving the formula for the mean first- 
passage time in the following text. 

Formula for the mean first-passage time. Having ob- 
tained the scaling of mean transmit time for old nodes, 
we now determine the mean first-passage time, with an 
aim to derive an exact solution. We represent the set of 
nodes in Ag as V,g, and denote the set of nodes created at 
generation g by Qg. Thus we have fig = fig U ^g-i. For 
the convenience of computation, we define the following 
quantities for m < g: 

-^m, total ~ 5^ -^i ' (^) 

and 

^m, total ^ ^ "^i- (10) 

Then, we have 

_ rpg . 7p9 / 1 1 \ 

^g, total ~ total ^ ^ g, total" V-^-^^ 

Next we will explicitly determine the quantity T^ To 
this end, we should firstly determine T^g toted- 
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We examine the mean time to absorption for the first 
several generations of Apollonian network. In the case of 
(7 = 1, by the very construction of Apollonian network, it 
follows that 



n = 
n - 
n = 
n = 



(l + Ti)/3- 
[l + TD/'d- 

(i + ri)/3- 

{l + Tl)/2>- 



(1 
(1 
(1 
(1 



n)/i- 



(1 + ^3^)73, 

(i + ri)A 

(l + Ti)/3. 



Thus, 



1, total 



= 4 + (Ti +Ti+Ti)^A + T,^,,,,, . (12) 
Similarly, for g = 2 case, it is easy to obtain 



20 



2, total 



4x3 + 2 Tq jQtjji + jQtjji , 



(13) 



Proceeding analogously, it is not difficult to derive that 



5, total 



= 4 X 3^ " 



t? — l, total 

93-2 7pS 
^ ^ 1, total 



2 2^g-2, total + ■ • ■ 
2^"'r^,total(14) 



and 



7^9+1 - 4 X 3f 

^ g+l, total — X O 



5, total 



2T 



p—l, total 



2^ ^ ^1, total + 2^ ^qj-q^jjI, (15) 



where 4 x 3^"^ and 4 x 3^ are actually the numbers of 
nodes generated at generations g and g + 1 , respectively. 
Equation (15) minus equation (14) times ^ and making 
use of the relation Tf^^ = | Tf , one gets 



- 5+1, total 



27 



T 



total 



X 39. 



(16) 



Using Tj (.Qjg^i equation (16) is solved inductively 



5, total 



32 (21 



39- 



(17) 



Substituting equation (17) for j^^j^j into equation (11) 
and using T^_i^t^,^i = § T^lltot^i, we have 

rpg _ 9 9_i , ^/^^y ,09-1 

^3, total g total ^ 9 



(18) 



Considering the initial condition T^^'^totai ^ equa- 
tion (18) is resolved by induction to yield 

(17x39 + 5^+1+32x99). (19) 



(/.total 



6 X 59 



V 



|~| Analytical Results 
# Numerical Results 



H 



N 



Fig. 6: Mean first-passage time {T)g versus network order A^g 
on a log-log scale. The solid line is a guide to the eye. 

Plugging the last expression into equation (7), we arrive 
at the accurate formula for the average of the mean time 
to absorption at the trap located at node 1 on the gth of 
Apollonian network: 



{T)g^ 



39(17 X 39 



32 X 99) 



6 X 59(2 X 39 + 1) 



(20) 



We have checked our analytic formula against numerical 
values quoted in Table 1. For the range of < g < 6, the 
values obtained from equation (20) completely agree with 
those numerical results on the basis of the direct calcula- 
tion through equation (7). This agreement serves as an 
independent test of our theoretical formula. 

From equation (1), we have g ~ logg (^ — l). Hence, 
for large network (i.e., Ng — > 00), we obtain 



- 1 



^(2-ln5/l„3)^ (21) 



where the exponent 2 — ~ 0.535 < 1. Thus, in the 
large limit of network order Ng, the MFPT increases al- 
gebraically with increasing order of the network. 

Generally, finite-size effect plays a crucial role in critical 
phenomena and dynamical processes on networks [45] . In 
order to study the finite-size effect in the scaling behavior 
of MFPT for the Apollonian network, we plot a figure to 
show the MFPT with order for small networks, see fig- 
ure 6. Surprisingly, it is observed from figure 6 that the 
finite-size effect has little influence on the scaling. 

Previous studies showed that for regular lattices [10], 
Sierpinski fractals [15,16], and T-graph [17] with large or- 
der TV, their MFPT (T) behaves as (T) - TV" with a > 1. 
Even in a complete graph of N nodes, Kn, the MFPT 
(T) grows linearly with TV as (T) = - 1. Actually, 
linear scaling of the MFPT with N is the best that has 
been reported so far. From equation (21), one can see 
that the MFPT (T) of the Apollonian network increases 
as a fractional power of network order iV, which implies 
that Apollonian network has a faster transmit time than 
any other analytically soluble media. In other words, the 
Apollonian network has the best structure for fast diffu- 
sion one can see heretofore. 
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We argue that the heterogeneous topology of the Apol- 
lonian network may be responsible for the high efficiency 
for transport. In the Apollonian network, there are some 
nodes with high degree, which are linked to one another 
and to most other nodes in the network. These 'large' 
nodes can be easily visited by a walker starting from an 
arbitrary location. Since 'large' nodes, including the trap 
node, are connected to one another, so the walker can find 
the trap in a very short time. 

Conclusions. — In summary, we have studied trap- 
ping problem on the Apollonian network exhibiting re- 
markable features (scale-free behavior and small-world ef- 
fects) of a variety of real networks. We have obtained an 
analytical closed-form solution for the MFPT for random 
walks with a trap located at a hub node, which is consis- 
tent with the numerical computation. The rigorous result 
indicates that the MFPT on the Apollonian network shows 
a very distinct behavior, which is compared with those re- 
sults previously reported for regular lattices, Sierpinski 
fractals, T-graph, and complete graph. We have shown 
that Apollonian network is the most efficient network for 
transport by diffusion, in contrast to other studied struc- 
ture. We hope that the current study may cast some light 
on trapping problem in some real networks, that present 
similar topologies as the Apollonian network. 
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